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Possible hadronization of supercooled QGP, created in heavy ion collisions at RHIC 
and SPS, is discussed within a Bjorken hydrodynamic model. Such a hadronization 
^ is expected to be a very fast shock-like process, what, if hadronization coincides or 

shortly followed by freeze out, could explain a part of the HBT puzzle, i.e. the flash- 
like particle emission {R ou t/Rside ~ 1)- HBT data also show that the expansion time 
before freeze out is very short (~ 6 — 10 fm/c). In this work we discuss question of 
supercooled QGP and the timescale of the reaction. 
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1 Introduction 

Two-particle interferometry has become a powerful tool for studying the size 
and duration of particle production from elementary collisions (e + e~, pp and 
pp) to heavy ions like Au + Au at RHIC or Pb + Pb at SPS [1,2]. For the 
case of nuclear collisions, the interest mainly focuses on the possible transient 
formation of a deconfined state of matter. This could affect the size of the 
region from where the hadrons (mostly pions) are emitted as well as the time 
for particle production. 

Comparing recent data [3] from RHIC with SPS data one finds a "puzzle" [4]: 
all the HBT radii are pretty similar although the center of mass energy is 
changed by an order of magnitude. Discussions at "Quark Matter 2002" [5] 
lead to the conclusion that the duration of particle emission, as well as the life- 
time of the system before freeze out, appear to be shorter than the predictions 
of most of the model at the physics market. 

It was demonstrated that a strong first-order QCD phase transition within 
continuous hydrodynamical expansion would lead to long lifetimes of the par- 
ticle source [2,6,7]Q, which would manifest itself as a large R ou t/R S ide ratio. 
Now this type of hadronization is excluded by experimental data. 

An alternative possibility, discussed in Refs. [9,10,11,12,13], is the hadroniza- 
tion from the supercooled QGP. This is expected to be a very fast shock-like 
process. If the hadronization from supercooled QGP coincides with freeze out, 
like it was assumed in Ref. [12], then this could explain a part of the HBT 
puzzle, i.e. the flash-like particle emission (R out / R S id e ~ 1)- In this work we are 
asking the following question - can the hadronization from supercooled QGP 
explain also the another part of the HBT puzzle, i.e. a very short (~ 6 [14] 
— 10[5,15] fm/c) expansion time before freeze out? 



2 Shock hadronization of a sQGP 

Relativistic shock phenomena were widely discussed with respect to their con- 
nection to high-energy heavy ion collisions (see, for example, [16]). In thermal 
equilibrium by admitting the existence of the sQGP and the superheated 
hadronic matter (HM) we have essentially richer picture of discontinuity-like 
transitions than in standard compression and rarefaction shocks. The system 



If we use some microscopic model for hadronization, for example nucleation 
of relativistic first-order phase transition [8], the lifetime is even longer - it was 
estimated to be about 50 — 100 fm/c in Ref. [8]. 
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evolution in relativistic hydrodynamics is governed by the energy-momentum 
tensor T^ v = {e-\-p)u^u v —pg^ v and conserved charge currents (in our applica- 
tions to heavy ion collisions we consider only the baryonic current n-u M ). They 
consist of local thermodynamical fluid quantities (the energy density e, pres- 
sure p, baryonic density n) and the collective four- velocity = — v 2 (l, v). 
Continuous flows are the solutions of the hydrodynamical equations: 



with specified initial and boundary conditions. These equations are nothing 
more than the differential form of the energy-momentum and baryonic number 
conservation laws. Along with these continuous flows, the conservation laws 
can also be realized in the form of discontinuous hydrodynamical flows which 
are called shock waves and satisfy the following equations: 



where da^ is the unit 4-vector normal to the discontinuity hypersurface. In 
eq. (2) the zero index corresponds to the initial state ahead of the shock 
front and quantities without an index are the final state values behind it. A 
general derivation of the shock equations (valid for both space-like and time- 
like normal vectors da^) was given in Ref. [17]. 

The important constraint on the transitions (2) (thermodynamical stability 
condition) is the requirement of non- decreasing entropy (s is the entropy 
density) : 

su^da^ > SoU^da^ . (3) 

To simplify our consideration and make our arguments more transparent we 
consider only one-dimensional hydrodynamical motion. To study the shock 
transitions at the surface with space-like (s.l.) normal vector (we call them s.l. 
shocks) one can always choose the Lorentz frame where the shock front is at 
rest. Then da^ = (0, 1) at the surface of shock discontinuity, and eq. (2) in 
this (standard) case becomes: 



(1) 



T^da v = T^dov , n u%da^ = nu^da. 



(2) 




ii 



n Q ul = nu 



i 



(4) 



Solving eq. (4) one obtains 




(p-Po)(t + Po) 
(e - e )(e + p) 



(p-Po)(to + p) 
(e - e )(e + p D ) 



(5) 
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Fig. 1. Possible final states in the (energy density-pressure)-plane for shock transi- 
tions from the initial state (e ,p ). I and IV are the physical regions for s.l. shocks, 
III and VI for t.l. shocks. II and V are unphysical regions for both types of shocks. 
Note, that only states with p < e are possible for any physical Equation of State in 
the relativistic theory. 

and the well known Taub adiabat (TA) [18] 

n 2 X 2 - n 2 X 2 - (p - Po ) (X + X ) = , (6) 



where X = (e + p)/n 2 !^\ 

For discontinuities on a hypersurface with a time-like (t.l.) normal vector da^ 
(we call them t.l. shocks) one can always choose another convenient Lorentz 
frame ("simultaneous system") where da' 1 = (1,0). Equation (2) is then 



2 It has been shown in a series of works [19], that freeze out through the space-like 
hypersurface leads to nonequilibrium post FO distribution. 



4 



Solving eq. (7) we find 

~2 (e- e )(e + p) 2 (e-e )(e + p ) 

v o = 7 77 7? v = 7 77 7 ? l°J 

{p-Po){e + Po) {p-Po){e +P) 

where we use the " ~ " sign to distinguish the t.l. shock case (8) from the 
standard s.l. shocks of (5). Another relation contains only the thermodynam- 
ical variables. It appears to be identical to the TA of eq. (6). Eqs. (8) and (5) 
are connected to each other by simple relations [11]: 

% = \ = \ ■ (9) 

These relations show that only one kind of transition can be realized for a 
given initial state and final state. The physical regions [0, 1) for v 2 ,v 2 (5) and 
for v 2 Q) v 2 (8) can be easily found in (e-p)-plane [11]. For a given initial state 
(e ,p ) they are shown in Fig. 1. For supercooled initial QGP states the TA 
no longer passes through the point (e ,p ) and new possibilities of t.l. shock 
hadronization transitions to regions III and VI in Fig. 1 appear. 



3 Hadronization of the sQGP within Bjorken hydrodynamics 



For a study of the expanding QGP we have chosen a framework of the one 
dimensional Bjorken model [20] (actually our principal results will not change 
if we use 3D Bjorken model). Within the Bjorken model all the thermodynam- 
ical quantities are constant along constant proper time curves, r = Vt 2 -z 2 = 
const. The important result of Bjorken hydrodynamics (which assumes a per- 
fect fluid) is that the evolution of the entropy density, is independent of the 
Equation of State (EoS), namely 

s(r) = . (10) 



In Bjorken model the natural choice of the freeze out hypersurface is r = const 
hypersurface, where normal vector is parallel to the Bjorken flow velocity, 
v = z/t. Thus, do v = (1, 0) in the rest frames of each fluid element. This leads 
to the simple solution of the t.l. shock equations (7): 

v 2 = vl = 0, e = e ,n = n , p ^ p Q . (11) 



The entropy condition (3) is reduced to 

s > s . (12) 
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Fig. 2. Different ways for a system to go from Q state (sq) to H state (sh) are pre- 
sented on {s, t} plane. Subplot A shows continuous expansion, which takes time th, 
eq. (13). Subplot B presents flash-like particle emission, i.e. simultaneous hadroniza- 
tion and freeze out; which takes time eq. (14). Subplot C shows several possi- 
bilities according to scenario 2 with shock-like hadronization into superheated HM. 
Time rff (16) can be smaller or larger than depending on details of the EoS, 
but always larger than th- 



Now let us try to answer the main question of this work - can QGP expansion 
with t.l. shock hadronization of supercooled state be faster than the hadroniza- 
tion through the mixed phase? The initial state is given at the proper time 
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Tina = tq, when the local thermal equilibrium is achieved in the QGP state 
Q = (cq,Pq, sq). The final equilibrium hadron state is also fixed, by experi- 
ment or otherwise, as H = (e H ,p H , s#). For the continuous expansion given 
by eq. (10) the proper time for the Q — > H transition is (see Fig. 2 - subplot 
A): ' 

TH = ^ . (13) 

sh 



If our system enters the sQGP phase and the particle emission is flash-like, i.e. 
the system hadronizes and freezes out at the same time, then eq. (10) is also 
valid all the time with final t.l. shock transition to the same H state. We call 
this as a scenario number one (see Fig. 2 - subplot A). Our system should go 
into supercooled phase to the point where = e#, nffl = nn, as it is required 
by eq. (11). At this point our sQGP has entropy density s£\ It's value depend 
on the EoS, but t.l. shock transition is only possible if < sh according 
to eq. (12). Thus, for the proper time of Q — > H transition according to first 
scenario we have: 

(i) S Q T Q ^ ( i A \ 
Th = — (TP > t h- (14) 

So 



We can also study a scenario number two when our system supercools to the 
state {t 2 \p 2 \ s^) : then hadronizes to a superheated HM state (e^ ,p^ 2 \ s^), 
and then this HM state expands to the same freeze out state H = (€h,Ph, sh)- 
(see Fig. 2 - subplot C). At the point of the shock transition one has: 

= (15) 

So 



Then we have a t.l. shock transition satisfying eq. (11), and following the HM 
branch of the hydrodynamical expansion we find: 

(2) * (2) sqTq 

th = = > t h, 16 

sh s h si' 



since > due to non-decreasing entropy condition (12). In this second 
scenario the value of the of entropy density of sQGP can be both smaller 
and larger than HM final value sh- Depending on details of the EoS the proper 
time (16) of the Q — > H transition can also be smaller as well as larger 
than (14), but always larger than t h . 
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4 Conclusions 



The conclusion of our analysis seems to be a rather general one: the system's 
evolution through a supercooled phase and time-like shock hadronization can 
not be shorter than a continuous expansion within the perfect fluid hydro- 
dynamics independently of the details of EoS and the parameter values of 
the initial, Q, and final, H, states. Although in such a way we may achieve 
flash-like particle emission, supported by the HBT data, the expansion time 
becomes longer, making it harder to reproduce the experimental HBT radii. 

So, how can we achieve shorter freeze out time than the minimal one com- 
ing from (very fast) Bjorken expansion via thermal and phase equilibrium? 
Any delay in the phase equilibration (see assignment 9 in Ref. [21]) or/and 
any dissipative process in our system lead to the entropy production, what in- 
creases the time needed to reduce entropy density to s Q < sh (for the flash-like 
particle emission). 

The system may, nevertheless, freeze out and hadronize into a non-equilibrated 
hadron gas well before th- This is possible e.g. through a dominantly s.l. 
hypersurface with non-decreasing entropy condition, eq. (12) [19], at earlier 
times from a slightly supercooled QGP. On the other hand a dominantly s.l. 
hypersurface gives a finite duration of the particle emission, making it harder 
to reproduce experimental R ou t/R S ide ratio. 

The construction of a full reaction model, which simultaneously describes data 
on two particle interferometry, hadron spectra and hadron abundances is a 
formidable task which is still ahead of us. 



Acknowledgements 

We have benefitted from discussions with I. Mishustin, L. Satarov and Yu.M. 
Sinyukov. The financial support from the Humboldt Foundation and INTAS 
grant 00-00366 are acknowledged. 



References 

[1] E. Shuryak, Phys. Lett. B44 387 (1973); M. Gyulassy, S. K. Kauffmann and 
L. W. Wilson, Phys. Rev. C 20 2267 (1979); Yu.M. Sinyukov, Nucl. Phys. A498 
151 (1989), A556 589 (1994); U. W. Heinz and B. V. Jacak, Ann. Rev. Nucl. 
Part. Sci. 49 529 (1999); R. M. Weiner, Phys. Rept. 327 249 (2000). 



8 



[2] S. Pratt, Phys. Rev. Lett. 53 1219 (1984); Phys. Rev. D 33 1314 (1986); 
G. Bertsch, M. Gong, M. Tohyama, Phys. Rev. C 37 1896 (1988); G. Bertsch, 
Nucl. Phys. A498 173 (1989). 

[3] RHIC data: C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 87 082301 

(2001) ; K. Adcox et al. [PHENIX Collaboration], Phys. Rev. Lett. 88 192302 

(2002) ; SPS data: S. V. Afana sev et al. [NA49 C ollaboration], Nucl. Phys. A698 
104 (2002); H. Appelshauser, |hep-ph/020415S| . 

[4] M. Gyulassy, jnucl-th/010607^ . 

[5] S. Pratt , talk at the Quark Matter 2002, July, Nantes, France. 
[6] D. Rischke, M. Gyulassy, Nucl. Phys. A608 479 (1996). 



[7] D. Teaney, J. Lauret,, E.V. Shuryak, Phys. Rev. Lett. 86 4783 (2001) and |nucl- 
| th/0110037| . 

[8] L.P. Csernai, J.I. Kapusta, Phys. Rev. Lett. 69 737 (1992); Phys. Rev. D 46 
1379 (1992). 

[9] T. Csorgo, L.P. Csernai, Phys. Lett. B333 494 (1994). 



[10 
[11 

[12 

[13; 

[14 
[15 

[ie; 

[17; 
[is; 

[19 

[20; 

[21 



L.P. Csernai, I.N. Mishustin, Phys. Rev. Lett. 74 5005 (1995). 

M.I. Gorenstein, H.G. Miller, R.M. Quick, R.A. Ritchie, Phys. Lett. B340 109 
(1994). 

A. Keranen, J. Mannien, L.P. Csernai, V. Magas, nucl-th/0205019| . 
L.L. Jenkovszky, B. Kampfer, V.M. Sysoev, Z. Phys. C 48 147 (1990). 
T. Csorgo, A. Ster, |nucl-th/02070T6 . 



S. Panitkin, talk at the XXXII International Sumposium for Multiparticle 
Dynamics, 7-13 Sep 2002, Alushta, Ukraine. 

R.Stock, Phys. Rep. 135 259 (1986); H.Stocker and W.Greiner, Phys. Rep. 137 
277 (1986); KA. Bugaev, M.I. Gorenstein, B. Kampfer and V.I. Zhdanov, Phys. 
Rev. D 40 2903 (1989). 

L. Csernai, Zh. Eksp. Teor. Fiz. (Russ.) 92 379 (1987); Sov. Phys. JETP 65 
216 (1987). 

A.H. Taub, Phys. Rev. 74 328 (1948). 

Cs. Anderlik et al, Phys. Rev. C 59 388 (1999); Phys. Rev. C 59 3309 (1999); 
Phys. Lett. B459 33 (1999); V.K. Magas et al, Heavy Ion Phys. 9 193 (1999); 
Nucl. Phys. A661 596 (1999). 

J.D. Bjorken, Phys. Rev. D 27 140 (1983). 

L.P. Csernai, Introduction to Relativistic Heavy Ion Physics, Wiley, 1994. 



9 



